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1 INTRODUCTION 


All flutter suppression systems require sensors to detect the movement of 
the lifting surface and to activate a control surface according to a synthe- 
sized control law. Most of the work performed to date (refs. 1 through 
5) relate to the development of control laws based on predetermined loca- 
tions of sensors and control surfaces. These locations of sensors and control 
surfaces are determined either arbitrarily, or by means of a trial and error 
procedure (ref. 5) 

The aerodynamic energy concept (ref. 6) indicates that the sensors 
should be located within the activated strip. Furthermore, the best chord- 
wise location of a sensor activating a T.E. control surface, is around the 
65-percent chord location (ref. 7). The best chordwise location 
for a sensor activating a L.E. surface is shown to lie upstream of the wing 
(around 20-percent upstream of the leading edge), or alternatively, two 
sensors located along the same chord should be used. 


Plan View of a Wing with Active Control Surfaces 
and Recommended Sensor Locations 




The present work describes a method which enables one to determine the 
best spanwise placement of an activated control surface without resorting to 
any specific control law. The method is based on the aerodynamic energy 
concept whereby the activated control surface is placed at the location 
where most energy is fed into the unstable structure. 

2 APPROACH 

Let the pressure p( x, y ) be given by eq.(l) and let the displacement z(x, y ) 
be given by eq.(2), where the q ' s denote the generalized coordinates of the 
system. The generalized aerodynamic forces per unit area are given by 
eq.(3), or in a more condensed form by eq.(4). 
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The work per unit area W x>y done by the system on its surroundings 
per cycle of oscillation is given by (ref. 6) eq.(5), where eq.(6) provides 
a definition for some of the parameters in eq.(5). If W Xi y > 0, energy is 
dissipated by the system. If W x>y < 0, energy is fed into the system. Eq.(8) 
is obtained by integrating eq.(5) along the chord, and it yields the work 
W y per unit span done by the system on its surroundings per cycle of 
oscillation. The total work W is obtained by integrating eq.(8) along the 
span, as shown in eqs.(10),(ll). 
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Define the specific energy ratio by W A , as given by eq.(13), and note 
that its integral along the span must have a unit absolute value (see eq.(14)). 


W A = 


Wy_ 
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It is argued that the best spanwise placement of an active control surface 
for flutter suppression is around the location where W A is negative and 
assumes the largest numerical values. 


3 DETAILED PROCEDURE 

1. Determine the flutter dynamic pressure Qp of the system 

2. Increase Qd , so as to lie within the unstable region, and obtain the 
eigenvector {go} of the unstable mode. The amount by which Qd is 
increased is immaterial since only energy ratios are used. 

3. Compute W (which must be negative), and W A for the different span- 
wise locations. 

4. Plot the specific energy ratio W A versus the span and determine the 
strip where W A assumes the largest negative value. This strip absorbs 
most energy per unit span and therefore would present the location 
where an active control surface would be most effective in suppressing 
flutter. 



4 NUMERICAL EXAMPLES 

Three numerical examples are used.These examples relate to 

1. DAST-ARW2,12 modes (2 rigid body modes), see ref. 8. 

2. OBLIQUE WING, 20 modes(5 rigid body modes), see ref. 9. 

3. MODIFIED OBLIQUE WING (with one torsional modal frequency 
reduced so as to cause wing flutter) 

The relevant aerodynamic matrices were computed using Langley’s dou- 
blet lattice ISAC program.The results obtained are shown in the following 
figures. 


Oblique Wing Model Planform View with Wing 
Skewed 65° 
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5 Results for the DAST-ARW2 Mathemat- 
ical Model 

The geometrical layout of the DAST-ARW2 model is shown below. It can 
be seen that out of the 24 strips allowed for the model, 17 strips lie along the 
wing, 4 strips lie along vertical tail and 3 strips lie along the horizontal tail. 
This model yields a flutter dynamic pressure Qf = 490 psf (at M=0.85) 
and a flutter frequency u>p = Wlrad/s. The unstable eigenvector was 
computed for Q = 550 psf and the matrices [A] and [A] y were computed 
for all the 24 strips at the reduced frequency k = 0.132 associated with the 
unstable mode. 

ARW2 Geometric Layout, Together with Doublet 
Lattice Paneling and Strip Number Allocations 








The figure below shows a plot of the specific energy ratios W a for only 
those 17 strips that lie along the wing. The values of W a are negligible 
for all the other strips and therefore will not be shown herein. As can 
be seen, the specific energy ratio is negative for all wing strips, except for 
the root strip (strip 5) where W a is very small and positive. The largest 
negative numerical value of W a relates to strip 18, which coincides with the 
inboard portion of the aileron. Following the method described herein, this 
is the location around which the aileron should be placed for best effects 
regarding flutter suppression ( i.e. around the 80-percent span location). 


DAST - ARW2 Model - Variation of Specific Energy 
Ratio W A with Strip Locations Along the Wing 



When W A >0, strip is 
dissipating energy 
When W A < 0, strip Is 
absorbing energy 
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6 Results for the Oblique Wing 

The geometrical layout of the wing in a 65-degrees skew position (with right 
wing forward) is shown. Note that the wing has again 17 strips along its 
span, with strip 6 at the tip of the left wing and strip 22 at the tip of the 
right wing. 


Oblique Wing Model - Geometric Layout Together 
with Doublet Lattice Paneling and Strip Numbers 

Planeform 65° Skew 
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The horizontal tail, the horizontal projection of the fuselage, and the 
side view of the vertical tail are both shown in the figure below. There are 
4 strips on each of the horizontal tails (left and right surfaces), 2 strips on 
the horizontal projection of the fuselage, and 4 strips on the vertical tail. 

A flutter computation at Mach 0.95 shows that a mild, 78 rad/s, vertical 
tail flutter instability develops around Qf = 780 psf. 


Oblique Wing Model - Geometric Layout Together 
with Doublet Lattice Paneling and Strip Numbers 



Vertical tail 
sideview 
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The specific energy ratio distribution was computed for Q=1600 psf. 
It can be seen that most of the energy input into this fluttering system 
takes place through the vertical tail and the tips of the horizontal tail. The 
inboard parts of the horizontal tail, the horizontal fuselage, and practically 
all of the wing, all dissipate energy and thus contribute to the mildness of 
the flutter obtained. The following figures indicate that for the suppression 
of this flutter mode, the active control surface should be placed around 
the center of strip 3 of the vertical tail (i.e. around its 60-percent span 
location). 

Oblique Wing Model - Variation of Specific Energy 
Ratio W A with Strip Locations 

Wing 



Span Y along wing, in. 
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Oblique Wing Model- Variation of Specific Energy 
Ratio with Strip Locations 

Horizontal Tail and Horizontal Fuselage 
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When W A >0, strip is 
dissipating energy 
When W A < 0, strip is 
absorbing energy 


Oblique Wing Model - Variation of Specific Energy 
Ratio W^ with Strip Locations 

Vertical Tail 



When W A >0, strip is 
dissipating energy 
When W A <0, strip is 
absorbing energy 



7 Results for the Modified Oblique Wing 

The oblique wing is unique from the point of view of assymetry. The best 
placement of an activated strip along its span could present an interesting 
challenge. The mathematical model of the wing was therefore modified so as 
to ‘force’ the wing to flutter. This was done by lowering one of the torsional 
frequencies of the wing from 45 HZ to 12 HZ. Therefore, the following results 
do not relate to the actual wing but to a synthetically modified wing. This 
modified model yields, in addition to the already seen vertical tail flutter, a 
wing flutter mode with Qf — 1050ps/ and Wf = 70 rad/ s. The following 
results relate to this wing flutter mode, with specific energy ratios computed 
at Q=1300 psf. 


Modified Oblique Wing Modef - Variation of 
Specific Energy Ratio W A with Strip Locations 
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Modified Oblique Wing Model- Variation of 
Specific Energy Ratio W A with Strip Locations 

Vertical Tail 



When W A >0, strip is 
dissipating energy 
When W A < 0, strip is 
absorbing energy 


Modified Oblique Wing Model- Variation of 
Specific Energy Ratio W A with Strip Locations 
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As seen from the above 3 figures, most of the energy input into the system 
takes place through the left wing. The vertical tail essentially dissipates 
energy, and the horizontal tail absorbs energy through its left surface and 
dissipates around the same amount of energy through its right surface. 
The largest negative numerical value for W a is obtained in strip 6 which 
represents the tip of the left (aft) wing. Hence, for the suppression of this 
flutter mode, the activated control surface should be placed as close to the 
left tip of the wing as is structurally possible. Strips 9 and 10 can also 
form a reasonable alternative to the above extreme tip placement of the 
activated control surface, i.e., around the 65-percent span location of the left 
wing. 
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